Electronic Properties of Disordered Two-Dimensional Carbon 
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Two-dimensional carbon, or graphene, is a semi-metal that presents unusual low-energy electronic 
excitations described in terms of Dirac fermions. We analyze in a self-consistent way the effects of 
localized (impurities or vacancies) and extended (edges or grain boundaries) defects on the elec- 
tronic and transport properties of graphene. On the one hand, point defects induce a finite elastic 
lifetime at low energies with the enhancement of the electronic density of states close to the Fermi 
level. Localized disorder leads to a universal, disorder independent, electrical conductivity at low 
temperatures, of the order of the quantum of conductance. The static conductivity increases with 
temperature and shows oscillations in the presence of a magnetic field. The graphene magnetic sus- 
ceptibility is temperature dependent (unlike an ordinary metal) and also increases with the amount 
of defects. Optical transport properties are also calculated in detail. On the other hand, extended 
defects induce localized states near the Fermi level. In the absence of electron-hole symmetry, these 
states lead to a transfer of charge between the defects and the bulk, the phenomenon we call self- 
doping. The role of electron-electron interactions in controlling self-doping is also analyzed. We 
also discuss the integer and fractional quantum Hall effect in graphene, the role played by the edge 
states induced by a magnetic field, and their relation to the almost field independent surface states 
induced at boundaries. The possibility of magnetism in graphene, in the presence of short-range 
electron-electron interactions and disorder is also analyzed. 
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I. INTRODUCTION 



Carbon is a life sustaining element that, due to the 
versatility of its bonding, is present in nature in many 
allotropic forms. Besides being an element that is 
fundamental for life on the planet, it has been ex- 
plored recently for basic science and technology in the 
form of three-dimensional graphite, 1 one-dimensional 
nanotubes^ zero-dimensional fullerenes^ and more re- 
cently in the form of two-dimensional Carbon, also known 
as graphene. Experiments in graphene-based devices 
have shown that it is possible to control their elec- 
trical properties by the application of external gate 
volt age ^*S*LSi2ii2iii opening doors for carbon-based 
nano-electronics. In addition, the interplay between dis- 
order and magnetic field effects leads to an unusual 
quantum Hall effect predicted theoreticallyi 2 **^*^ and 
measured experimentally^*^. These systems can be 
switched from n-type to p-type carriers and show entirely 
new electronic properties. We show that their physical 
properties can be ascribed to their low dimensionality, 
and the phenomenon of self-doping, that is, the change in 
the bulk electronic density due to the breaking of particle- 
hole symmetry, and the unavoidable presence of struc- 
tural defects. Our theory not only provides a descrip- 
tion of the recent experimental data, but also makes new 
predictions that can be checked experimentally. Our re- 
sults have also direct implication in the physics of Carbon 
based materials such as graphite, fullerenes, and carbon 
nanotubes. 

Graphene is the building block for many forms of Car- 



bon allotropes. Its structure consists of a Carbon hon- 
eycomb lattice made out of hexagons (see Fig. The 
hexagons can be thought of Benzene rings from which 
the Hydrogen atoms were extracted. Graphite is ob- 
tained by the stacking of graphene layers that is sta- 
bilized by weak van der Waals interactions^ Carbon 
nanotubes are synthesized by graphene wrapping. De- 
pending on the direction in which graphene is wrapped, 
one can obtain either metallic or semiconducting elec- 
trical properties. Fullerenes can also be obtained from 
graphene by modifying the hexagons into pentagons and 
heptagons in a systematic way. Even diamond can be ob- 
tained from graphene under extreme pressure and tem- 
peratures by transforming the two-dimensional sp 2 bonds 
into three-dimensional sp 3 ones. Therefore, there has 
been enormous interest over the years in understanding 
the physical properties of graphene in detail. Neverthe- 
less, only recently, with the advances in material growth 
and control, that one has been able to study truly two- 
dimensional Carbon physics. 

One of the most striking features of the electronic 
structure of perfect graphene planes is the linear rela- 
tionship between the electronic energy, Ep. , with the two- 
dimensional momentum, k = (k x ,k y ), that is: E^ = 
Vp\k\, where vf is the Dirac-Fermi velocity. This sin- 
gular dispersion relation is a direct consequence of the 
honeycomb lattice structure that can be seen as two in- 
terpenetrating triangular sub-lattices. In ordinary metals 
and semiconductors the electronic energy and momen- 
tum are related quadratically via the so-called effective 
mass, m* , (Ek = ti 2 k 2 / (2m*)), that controls much of 
their physical properties. Because of the linear dispersion 
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relation, the effective mass in graphene is zero, leading 
to a unusual electrodynamics. In fact, graphene can be 
described mathematically by the two-dimensional Dirac 
equation, whose elementary excitations are particles and 
holes (or anti-particles), in close analogy with systems in 
particle physics. In a perfect graphene sheet the chemi- 
cal potential, /i, crosses the Dirac point and, because of 
the dimensionality, the electronic density of states van- 
ishes at the Fermi energy. The vanishing of the effective 
mass or density of states has profound consequences. It 
has been shown, for instance, that the Coulomb interac- 
tion, unlike in an ordinary metal, remains unscreened^ 
and gives rise to an inverse quasi-particle lifetime that in- 
creases linearly with energy or temperature^, in contrast 
with the usual metallic Fermi liquid paradigm, where the 
inverse lifetime increases quadratically with energy. 

The fact that graphene is a two-dimensional system 
has also serious consequences in terms of the positional 
order of the Carbon atoms. Long-range Carbon order in 
graphene is only really possible at zero temperature be- 
cause thermal fluctuations can destroy long-range order 
in two-dimensions (the so-called, Hohenberg-Mermin- 
Wagner theorem^). At a finite temperature T, topo- 
logical defects such as dislocations are always present. 
Furthermore, because of the particular structure of the 
honeycomb lattice, the dynamics of lattice defects in 
graphene planes belong to the generic class of kineti- 
cally constrained models22i2i, where defects are never 
completely annealed since their number decreases only as 
a logarithmic function of the annealing time 2 ^. Indeed, 
defects are ubiquitous in carbon allotropes with sp 2 co- 
ordination and have been observed in these systems^. 
As a consequence of the presence of topological defects, 
the electronic properties discussed previously, are signifi- 
cantly modified leading to qualitatively new physics. As 
we show below, extended defects can lead to the phe- 
nomenon of self-doping with the formation of electron or 
hole pockets close to the Dirac points. We show, how- 
ever, that the presence of such defects can still lead to 
long electronic mean free paths. We present next an anal- 
ysis of the physical properties of graphene as a function of 
the density of defects, at zero and finite temperature, fre- 
quency, and magnetic field. The defects analyzed here, 
like boundaries (edges), dislocations, vacancies, can be 
considered strong distortions of the perfect system. In 
this respect, our work complements the studies of de- 
fects and interactions in systems described by the two- 
dimensional Dirac equation 23 . 

The role of disorder on the electronic properties of 
coupled graphene planes shows also its importance on 
the unexpected appearance of ferromagnetism in proton 
irradiated graphito 24 i 25 i 26 i 27 i 28 i 29 . In a recent publica- 
tion, the role of the exchange mechanism on a disordered 
graphene plane was addressed 3 ^. It was found that disor- 
der can stabilizes a ferromagnetic phase in the presence 
of long-range Coulomb interactions. On the other hand, 
the effect of disorder on the density of states of a single 
graphene plane amounts to the creation of a finite density 



of states at zero energy. Therefore, a certain amount of 
screening should be present and the question of whether 
the interplay of disorder and short-range Coulomb inter- 
action may stabilize a ferromagnetic ground state has to 
be addressed as well. 

Moreover, with the current experimental techniques, it 
is possible to study not only a single layer of graphene 
but also graphene multi-layers (bilayers, trilayers, etc). 
Recent experiments provide direct evidence that while 
the high-energy physics of graphene multi-layers (for en- 
ergies above around 100 meV from the Dirac point) is 
quite different from that of single layer graphene, the low- 
energy physics seems to be universal, two-dimensional, 
independent of the number of layers, and dominated by 
disorder^Siii. Hence, the work described here maybe 
fundamental for the understanding of this low-energy be- 
havior. There is still an interesting question whether this 
universal low-energy physics survives in bulk graphite. 

In this paper we present a comprehensive and 
unabridged study of the electronic properties of graphene 
in the presence of defects (localized and extended), and 
electron-electron interaction, as a function of tempera- 
ture, external frequency, gate voltage, and magnetic field. 
We study the electronic density of states, the electron 
spectral function, the frequency dependent conductivity, 
the magneto-transport, and the integer and fractional 
quantum Hall effect. We also discuss the possibility of 
a magnetic instability of graphene due to short-range 
electron-electron interactions and disorder (the problem 
of ferromagnetism in the presence of disorder and long- 
range Coulomb interactions was discussed in a previous 
publication 3 ^) . 

The paper is organized as follows: in Sec. Ill Al a for- 
mal solution for the single impurity and many impurities 
T— matrix calculation is given. Details of the position av- 
eraging procedure are given in Sec. [3 in connection with 
the same problem, but in a magnetic field. In Sec. Ill B| the 
problem of Dirac fermions in a disordered honeycomb lat- 
tice is studied within the full Born approximation (FBA) 
and the full self-consistent Born approximation (FSBA) 
for the density of states. Using the results of Sec. Ill Bl 
the spectral and transport properties of Dirac fermions 
are computed in Sec. IIIII In Sec. \lV\ we address the 
question of magnetism and the interplay between short- 
range electron-electron interactions and disorder. The 
density of states of Dirac fermions in a magnetic field 
perpendicular to a graphene plane is studied in Sec. Ivl 
and the magneto-transport properties of this system are 
computed both at zero and finite frequencies, using the 
FSBA. The quantization values for the integer quantum 
Hall effect and for Jain's sequence of the fractional quan- 
tum Hall effect are discussed. Finally, Sec. IVIII contains 
our conclusions. We have also added appendices with the 
details of the calculations. 
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II. IMPURITIES AND VACANCIES. 

The honeycomb lattice can be described in terms of 
two triangular sub-lattices, A and B (see Fig. 0. The 
unit vectors of the underlying triangular sub-lattice are 



a a = -(3,V3,0), 
a 2 - !(3,-\/3,0), 



(2.1) 



where a is the lattice spacing (we use units such that 
Kb = 1 = h). The reciprocal lattice vectors are: 

6 1 = ^(1,73,0), b 2 = ^(1,-^3,0). (2.2) 
6a 6a 

The vectors connecting any A atom to its nearest neigh- 
bors are: 

Si = |(l,\/3,0), 

<* 2 = |(1,-V3,0), 
6 3 = a(l,0,0) (2.3) 
and the vectors connecting to next-nearest neighbors are: 



(2.4) 
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FIG. 1: (color on line) A honeycomb lattice with vacancies, 
showing its primitive vectors. 

In what follows we use a tight-binding description for 
the 7r-orbitals of Carbon with a Hamiltonian given by: 

Ht.b. = -t ]T (4 A- + h - c -) 

+ t' E (<.%,- A- + h " c -)' (2-5) 

where aj ^ (ai. CT ) creates (annihilates) and electron on site 
Ri with spin a (er =|, J,) on sub-lattice A and (&i !C7 ) 



creates (annihilates) and electron on site Hi with spin 
(7 (a =T)I) on sub-lattice B. t is the nearest neighbor 
hopping energy (t 2.7 eV), and t' is the next- 
nearest neighbor hopping energy (t'/t w 0.1). 
We notice en passant that in earlier studies of graphite 31 
it has been assumed that t' = 0. This assumption, how- 
ever, is not warranted since there is overlap between Car- 
bon 7r-orbitals in the same sub-lattice. In fact, we will 
show that t' plays an important role in graphene since it 
breaks the particle-hole symmetry and is responsible for 
various effects observed experimentally. 

Translational symmetry is broken by the presence of 
disorder. Localized defects such as vacancies and impu- 
rities are included in the tight-biding description by the 
addition of a local energy term: 



-ffimp. — ^ Vi 



(2.6) 



where Vi is a random potential at site R^. In momentum 
space we define: 

J _ 1 V ik-Rj t u\ _ 1 V" ifc-Rjfrt 

VN A k VN B k 

(2.7) 

where Na = Nb = N, and the non-interacting Hamilto- 
nian, Hi = if t .b, + i?im P ., reads: 



V l[ a h+q,a a k,a + 4+g,<A,<x] , (2. 

q,k,(7 



where 



k-(5< 



4>{k) = 

i=l 

6 

4>(k) = t'Y,e lkn 



(2.9) 



and Vq is the Fourier transform of the random potential 
due to impurities. Hamiltonian l|2.8|) is the starting point 
of our approach. 



A. The single impurity problem and the T-matrix 
approximation 

In the single impurity case one can write Vq = V/N 
where V is the strength of the impurity potential. In 
what follows we use standard finite temperature Green's 
function formalism 3 ^ 3 .. Because of the existence of two 
sub-lattices, the Green's function can be written as a 2 x 2 
matrix: 



G CT (fc,p,r) 



G A A,a{k,p,r) G A B,a(k,p,r) 

G B A,a{k,P,T) G B B,a(k,P,T) 



1(2.10) 
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where 



G A A, a(k,P,T) 

G AB , a{k,p,T) 
G B A : a(k,p,T) 



■(Ta k ,a(r)al t(T (Q)) , 
■(Ta fe)CT (r)&t CT (Q)) , 

■<T6 M (r)6t )jr (Q)), (2.11) 



where r is the "imaginary" time, and T is the time or- 
dering operator. 

For a single impurity the Green's function can be writ- 
ten as G(k, p, r) = Sk, P G(k,r), where2£: 

G(k, u n ) = G°(k, cu n ) + G°(k, u n )T imp .(u n )G (k, w n ) , 

(2.12) 

where u> n — 2irT(n+ 1/2) is the fermionic Matsubara fre- 
quency, G (fe, u) n ) is the propagator of the tight-binding 
Hamiltonian l|2.5|) and 

T imp .K) = ^[1 - VG ^)}- 1 , (2.13) 
is the single impurity T-matrix, where: 

G°K) = ^E G °( fc >^)- (2-14) 

The above result is exact for a single impurity. For a 
finite but small density, m = Ni/N, of impurities, the 
Green's function equation becomes: 

G(k, «„) = G°(fc, « n ) + G°{k, w n )T(w„)G(fc, c^„) , 

(2.15) 

which is valid up to first order in rij, that is, it takes only 
into account the multiple scattering of the electrons by a 
single impurity. Equation Q2.15JI can be solved as: 

G(k,u n ) = [[G?(k,u n )]- 1 -T(u n )]- 1 , (2.16) 

where 

T(Wn) = ^T imp .K) = Vm[l - VG ^)]- 1 . (2.17) 

For vacancies we take V — > oo and (|2.17(l reduces to: 

T(u n ) = -ml&fun)]- 1 . (2.18) 

It worth stressing that Eqs. 12.12fl and (I2.13[l although 
similar in form to Eqs. (|2.16(l and lj2.17|l have a very dif- 
ferent meaning. Whereas the first set applies to the sin- 
gle impurity problem, the latter set is the consequence 
of an assemble average over the impurity positions (see 
Sec. 0for details on the averaging procedure in the con- 
text of Landau levels) with a re-summation procedure, 
corresponding to the FBAS. 



B. The low-energy physics and the electronic 
density of states 

The results of the previous subsection are entirely gen- 
eral, in the sense that no approximation for the band 



structure was made. Consider, for simplicity, the tight- 
binding Hamiltonian Ij2.5|l in the case of t' = 0, that can 
be written, in momentum space, as: 



H, 



t.b. 



k,cr 



*L,<r>' 



<Xk) " 




<2k >( T 


>*(k) o 




&k, CT 



(2.19) 



which can be diagonalized and produces the spectrum: 

£±(k)=±|^(k)|, (2.20) 

where the plus (minus) sign is related with the upper 
(lower) band. It is easy to show that the spectrum van- 
ishes at the K point in the Brillouin zone with wave- 
vector, Q = (27r/(3V3a),27r/(3a)), and other five points 
in the Brillouin zone related by symmetry. In fact, it is 
easy to show that: 



0(Q + p) 



<MQ + p) 
I0(Q + p)I 



-tae l7T/3 {p y - ip x ) 



3 



e «5(Q+p) _ e , 



t/3 (Pv % Vx) | (2 22) 



where p (p -C Q) is measured relatively to the K point in 
Brillouin zone and we have def ined e l5{ ^ = <j>{k)/\4>{k)\, 
for latter use. Using (|2.21|) in l|2.20|l we find: 



E±(Q 



±\ta\p\ 



±v F \p\ 



(2.23) 



for the electron's dispersion. Eq. I|2.23|l is the dispersion 
of a relativistic particle with "light" velocity vf — 3ta/2, 
that is, a Dirac fermion. Hence, at low energies (energies 
much lower than the bandwidth), the effective descrip- 
tion of the tight-binding problem reduces the 6 points 
in the Brillouin zone to 2 Dirac cones, each one of them 
associated with a different sublattice. The low energy de- 
scription is valid as long as the characteristic momenta 
(energy) of the excitations is smaller than a cut-off, k c 
(D = vpkc), of the order of the inverse lattice spacing. 
In the spirit of a Debye model, where one conserves the 
total number of states in the Brillouin zone, we choose k c 
such that 7rfc^ = (2tt) 2 /A Ci where A c = 3\/3a 2 /2 is the 
area of the hexagonal unit cell. Hence, ea. (|2.23(l is valid 
for p <C k c and E <C D = k c VF- 

So far we have discussed the case of t' — 0. When t' ^ 
the problem can also be easily diagonalized and one finds 
that, close to the K point, the electron dispersion changes 
to: 



£±(Q + p) » -3t'±Hp| 



9t'a 2 



(2.24) 



showing that t' does not change the Dirac physics but 
introduces an asymmetry between the upper and lower 
bands, that is, it breaks the particle-hole symmetry. 
Hence, t' affects only the intermediate to high energy be- 
havior and preserves the low-energy physics. For many 
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of the properties discussed in this section t' does not play 
an important role and will be dropped. Nevertheless, we 
will see later that in the presence of extended defects t' 
plays an important role and has to be introduced in order 
to provide a consistent physical picture of graphcnc. 
For t' = we find: 



GAB(w n ,k) 

G B A(u n ,k) 



E 
E 



1/2 



ILJ„ 



E 

Gaa (w n , fc) 



j e «(fe) /2 
-«(*>) /2 



(2.25) 

(2.26) 

(2.27) 
(2.28) 



The expansion of the energy around the K point sim- 
plifies greatly the expressions in the calculation of the 
T-matrix, since they lead to simpler forms to Eqs. (|2.25|l - 
(|2.28(l . For the case of vacancies, Ea. l|2.18[l . it is easy to 
see that at low energies the T-matrix reads: 



T(oj n ) = -nj[G%(o; n )]" 1 1 
where 7 a 2x2 identity matrix, and 

G° AA (u> n ) = E 



(2.29) 



2N ^ ilo. 
j'=±l,k 



d 2 k 1 
(27r) 2 iuj„ — jvFk 

1 r k " dkk 



-E 



i 



Anrpvi 



iuj n — jvpk 
iu n ln(D 2 /cjl) , (2.30) 



where p = S/V is the graphene planar density (5 is 
the area of the graphene layer). After a Wick rotation 
(iuj n — ► lo + i0 + ) one finds: 



G%(w + i0+) = -F (u) - ZTrpoM . 



where, 



Po(w) 



IN 
D 2 



(2.31) 

(2.32) 
(2.33) 



where we have used that p = 1/A C = A: 2 /(47r), and hence 
Airpvp = D 2 . In the above equations we always assume 
\oj\ -C D. Notice that Po(oj) is simply the density of states 
of two-dimensional Dirac fermions. 



1. A single vacancy 

Assuming that a single unit cell has been diluted, we 
use Eqs. (|2.12|) and (|2.13|) to determine the correction to 



the Dirac fermion density of states. The actual density 
of states, p(ui), is given by: 

p(uj) = --ImG A A(u + iO + ) 



po(u) - 



2/N goM 
D 2 F 2 (u)+n 2 p 2 (u) 
2/N 



\u\\og z (D/\u\) 



(w -> 0) ,(2.34) 



indicating that the contribution of the vacancy to the 
density of states is singular in the low frequency regime. 
The contribution is negative because one has exactly one 
missing state associated with the vacancy. The elec- 
tronic wave function around a single impurity was com- 
puted in Ref . pjij . The result obtained here is identical to 
the one obtained in the dilution problem in Heisenberg 
antiferromagnetsi&Si. The reason for this coincidence 
is easy to understand: the low energy excitations of an 
antiferromagnet in the ordered Neel phase are antiferro- 
magnetic magnons with linear dispersion relation, that 
is, relativistic bosons with a "speed of light" given by 
the spin-wave velocity. Since we have been discussing a 
non-interacting problem, the statistics plays no role, and 
the effect of disorder is the same for relativistic bosons 
or fermions. 



2. The full Born approximation (FBA) 

The situation is clearly different if one has a finite 
density of vacancies. In this case we have to deal 
with Eqs. (|2.16|l and l|2.17|l corresponding to the FBA 
where all one-impurity scattering events have been con- 
sidered. As before, the density of states is given by 
p(ui) — —ImG aa{u + iQ + )/tt and it is possible, after 
some tedious algebra, to obtain an analytical expression 
for this quantity, given by: 



In 



D 2 



Po(uj) 2uj 
D 2 a(cu)^\b 2 (uj)+c 2 (Lu) 



ttD 2 



E 



arctan 



ab(u>) 
ab{u>) 



arctan 



a(w)D 



(2.35) 



with 



a{uj) = Fq(uj) + Tr 2 pl(u) , 
b(u>) — a(u>)w — niFo(u) , 
c(u>) = riiwpo(u}) , 



(2.36) 



where Fq(u) and pa(ui) are defined in (|2.32|l and l|2.33|l . 
respectively. A plot of Eq. (|2.35|l is given in Fig. for 
two values of the impurity concentration Uj. Once again, 
the low energy behavior of p{u>) is the same found in 
the context of diluted antiferromagnetsi2^& We remark 
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that the dilution procedure introduces a low energy scale 
proportional to Drii, as can be seen from panel (c) in 
Fig.0 
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FIG. 2: (color on line) Density of states obtained from the 
FBA. Panel (a) shows p(uj) over the entire band; panel (b) 
shows the low energy part, where it is seen that the peak in 
p{ui) has a higher value for lower m; panel (c) shows that the 
peak in p(ui) appears at an energy scale of the order of niD/4. 



3. The full self-consistent Born approximation (FSBA) 

The FBA does not take into account electronic scat- 
tering from multiple vacancies, it accounts only for mul- 
tiple scattering from a single one. In order to include 
some contributions from multiple site scattering another 
partial series summation can be performed by replacing 
the bare propagator in the expression of the T-matrix in 
(|2.18|> by full propagator, leading to the FSBA. Because 
the matrix elements of the scattering potential computed 
from two Bloch states \k) and \p) are assumed momen- 
tum independent, the self-energy for the electrons de- 
pends only on the frequency. The self-consistent problem 
requires, in general, a careful numerical solution but in 
this particular case it is possible to reduce the problem to 
a set of coupled algebraic equations. The self-consistent 
problem requires the solution of the equation: 



S(w„) = 



(2.37) 



where T,(uj n ) is the electron self-energy. One can show 
that the self-energy can be written as: 

E(^ + zO+)= — 1 , (2.38) 

b [u>) + nrp{uj) 

where F(u>) and p(u>) are determined by the following set 
of coupled algebraic equations: 



F(u) 



'" *(F,p,u;) + -^Lr(F,p,u;)(2.39) 



2a{uj)D 2 



i{uj)D 2 



= o TTFS *&P><>>) - ^LT(F,p, W )(,2.40) 



where we used the definitions i|2.36[l and also defined the 
functions ^f(F, p,ui) and T(F,p,w): 

'(aa(u)D + b(uj)) 2 + c 2 (uj)' 



y(F,p,u)J2^ 



(2.41) 



6 2 H +c 2 (w) 
T(F,p,u>) — — 2 arctan[6(o;)/c(u;)] 
^aarctan[a(w)D/c(w) - ab{u)/c(u>)} . (2.42) 



The solution of Eqs. I|2.39|l and (|2.4U|) describes the ef- 
fect of the vacancies on the density of states of the Dirac 
Fermions. p(u>) is the self-consistent density of states, 
and F(ui) corresponds to the real part of self-energy 
(in analogy with po(u>) and Fq(lo) defined in l|2.33|l and 
(|2.32|) ). In Fig. |3 we show the result of this procedure 
for various impurity concentrations. 
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FIG. 3: (color on line) Density of states obtained in FSBA. 
Panel (a) shows p(w); panel (b) shows F(u>); panel (c) shows 
p(ui) at low energies. 

The low-energy behavior of the density of states, show- 
ing a parabolic enhancement of p(uj), has also been found 
in the context of heavy-fermion superconductors^. An 
exact numerical calculation of the electronic density was 
carried out in Ref.j^j where it was found that besides 
the low energy dome-like shape of the p{uj) (as shown in 
Fig. |3J), a large peak appears very close to uj = 0. This 
peak is reminiscent of the single impurity result given in 
(|2.34l) . Hence, besides the peak, the FSBA gives a very 
good account of the density of states in this problem. 

Notice that the self-energy, S(w), in l|2.38p depends on 
rii in a non-trivial way, since both the self-consistent F(u) 
and p(oj) also depend on rii. The self-energy is depicted 
in Fig. 2] for various values of the dilution density rii. 



III. SPECTRAL AND TRANSPORT 
PROPERTIES 

The electronic spectral function is defined as: 
A(k,u) = --lmG(k,uj + iQ + ), 

7T 



(3.1) 



7 




FIG. 4: (color on line) Imaginary (left panel) and real (right 
panel) part of the self-energy obtained from the FSBA. Note 
that both quantities are divided by Dm. 



FIG. 5: (color on line) Imaginary part of the electron's self 
energy including both the effect of disordered and electron- 
electron interaction. The inset shows an intensity plot for the 
spectral function A(k,uj) for D = 8.248 eV, m = 0.0001. 



and can be interpreted as the probability density that 
an electron has momentum k and energy u>. For a non- 
interacting, non-disordered, problem, the spectral func- 
tion is simply a Dirac delta function at uj = E(k). In 
the presence of disorder and/or electron-electron interac- 
tions the spectral function is broadened and its sharpness 
determines whether the electronic system supports quasi- 
particles. The spectral function can be measured directly 
in angle resolved photoemission experiments (ARPES) 38 . 

In terms of the self-energy, E(k, oj), the spectral func- 
tion reads: 



A(k,u) 



1 



Im£(k,w) 



- E(k) - ReE(k,w)] 2 + [ImE(k,w)] 



43.2) 



In the case of graphcnc, there are two contributions to 
the self-energy, 



E(k,w) = Se.-e.(fe) + Edis.M , 



(3.3) 



where E c ._ G .(fc) is the self-energy correction due to the 
electron-electron interactions that was computed origi- 
nally in Ref. [13 : 



ImE r 



e (*) = — 
V ' 48 



v F \k\ , 



(3.4) 



where e is the electron charge, and eo the dielectric con- 
stant of graphene. The other contribution, E^is., is due 
to disorder and is given in Ij2.38|l . 

Notice that these two contributions to the self-energy 
have very different dependence with the energy: while the 
electron-electron self-energy decreases as the energy (mo- 
mentum) decreases, the self-energy due to disorder in- 
creases as the energy decreases. Hence, electron-electron 
interactions are dominant at high energies while disorder 
is dominant at low energies. This interplay between the 
two self-energies leads to the prediction that there will be 
a minimum in the self-energy for some energy where the 



electron-electron interaction becomes of the same order 
of the electron- vacancy interaction. In Fig. [S] we plot the 
self-energy as a function of energy for various impurity 
concentrations together with the spectral function (in- 
set). One can clearly observe the non-monotonic depen- 
dence of the self-energy with the energy. This behavior 
should be observable in ARPES experiments. 

Assuming an electric field applied in the ir-direction, 
the frequency dependent (real part) conductivity is cal- 
culated from the Kubo formula: 



1 f° 

o-(q,w) = - / 
w Jo 



d^ t ([4(q,t),J x (q,0)]) (3.5) 



where J x is the x-component of the current operator 
which, due to gauge invariance, has the formifi: 



J x = —ite u x ■ 8a\ a b i+ s^ - u x ■ 8b\ 



cr+8 L 



(3.6) 



(the notation i + S means Ri +6). In Fourier space, 
and after expanding the general expression around the 
K-point in the Brillouin zone, we obtain: 

K/3 bL^)- (3-7) 



J x = -iv F e 2 f (e~ 

k.a 



wr/3 t U 



Substitution of (|3.7|l into (|3.5(l shows that the problem 
depends on the Green's functions defined in 12.11fl . How- 
ever, due to the special form of Eq. I|2.22ll the conduc- 
tivity does not have contributions coming from products 
of Green's functions of the form GabGba- Taking into 
account the number of bands and the spin degeneracy, 
the Kubo formula for the real part of the conductivity at 
finite frequency and temperature has the form: 



<r(u,T) 



4v 2 p 2 f°° df 

*fc £lK« + „)-/(.)]x 



E 

k 



lmG A A(k, e + iO + )ImG AA (fc, e + w + i0 + ) , (3. 
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FIG. 6: (color on line) Left: Kernel of the integral for <r(0, T) 
as function of the energy. Product of the Fermi function 
derivative by K(e) at two different temperatures. 



FIG. 7: (color on line) Dependence of a(T, 0) on the temper- 
ature and on the impurity dilution n%. 



where /(e) = l/(e' e- ' 1 " !r + 1) is the Fermi-Dirac dis- 
tribution function. The integral over k in (|3.8() can be 
performed and find: 



a(u,T) 
where 



2tt 2 uj 



&[/( e + w)-/(e)]jr(w,<(fr9) 



K(u, e) = ImS(e + w) ImS(e) Q(w, e) , (3.10) 

(Q(u>, e) defined in Appendix IXjl . 

It is instructive to consider the zero-temperature, zero- 
frequency limit of the conductivity in Eq. (|3.8() (restoring 
h): 



2e 2 , 



[ImE(O)] 5 



[ImE(O)] 5 



2 e 2 

7T h 



The result l|3.11|l shows that as long as ImE(0) <C 
D, do has a universal value independent of the dilu- 
tion concentration, in agreement with earlier theoretical 
works 41,42 , and in agreement with the experimental data 
in graphena&. 

At finite temperatures the integral in (|3.9(l has to be 
evaluated numerically. Consider tr(0, T) whose behavior 
is determined by K(e) = K(0,e). The quantities K(e) 
and —f'(e)K(e) ( —/'(e) is the derivative of the Fermi 
function in order to e) are both represented in Fig. El The 
behavior of K(e) shows, "V"-like shape as the energy e 
is varied. As a consequence, er(0, 0) should present the 
same "V'Tike shape as the chemical potential fi moves 
around \jl = 0. Such behavior has indeed been observed 
in atomically thin carbon films^, where the density of 
electrons was controled by a gate potential. The temper- 
ature dependence of the <r(0, T), for fx = 0, is depicted 
in Fig. {7\ for different vacancy concentrations, and it is 
found to follow Sommerfeld asymptotic expansion, but 



the number of terms needed to fit the numerical curve 
grows very fast as the dilution is reduced. 

In Fig. |S] we plot the frequency dependence of a(u,T) 
obtained from numerical integration of (|3.9|) with the self- 
energy given in l|2.38|l . At low temperature, we see that 
a(u,T) develops a maximum around an energy value that 
is dependent on the number of impurities. In fact, if 
plot a(u>,T) as function of w/y/rii, the conductivity al- 
most shows scaling behavior for all impurity dilutions 
(see lower left panel). As the temperature increases, and 
if ct(0, T) is sufficiently large, the conductivity a(ui,T) 
acquires a Drude-like behavior (right panel). 
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FIG. 8: (color on line) Dependence of a(T,u>) (in units of 
e 2 /(nh)) , on the frequency w. 
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IV. MAGNETIC RESPONSE AND THE ROLE 
OF SHORT RANGE COULOMB INTERACTIONS 

The ferromagnetism measured in proton irradiated 
graphite opens the question whether the interplay of in- 
teractions and disorder can drive the system from a para- 
magnetic to a magnetic ground stated. We study the 
effect of disorder on the magnetic susceptibility in the 
presence of short-range interactions, and the resulting 
change in the tendency towards magnetic instabilities. 
The problem of magnetic instabilities due to long-range 
exchange interactions 43 in the presence of small density 
of carriers was discussed in great detail in ref. |3(j • We 
do not address here the effects associated to the interplay 
between the long-range Coulomb interaction and differ- 
ent types of lattice disorder— and the appearance of local 
moments close to defects^ 4 ^*^ 4 ^. 

The paramagnetic susceptibility of graphene is given 
by: 



X(T) = 



di 



d 1 



dh 



k n 



-4 



def'(e)p(e) , 



where 



(4.1) 



(4.2) 



is the magnetization, and p(e) is the electronic density 
of states which, in the presence of disorder, is given in 
(|2.40(l . Within the Stoner mechanism 4 -, ferromagnetism 
is possible if the local electron-electron interaction term 
(the so-called Hubbard term)^, U, is large than a critical 
value given by: 



(4.3) 



In the case of an antiferromagnetic instability the same 
criteria would lead to another critical value of U given 
by: 



1 2 f D 

— — = — / dep(e) arctan 



D 



ImS(e) 



(4.4) 



Notice that in the case of antiferromagnetism one finds 
that U% F ~ — rii) when ImE — > 0, in agreement 

with Hartree-Fock calculations^. 

In Fig.[!|]wc plot the magnetic susceptibility as function 
of T for different values of rij . The signature of the pres- 
ence of Dirac fermions comes from the linear dependence 
on T for T/D <C 1. Notice that, unlike the case of an 
ordinary metal that has a temperature dependent Pauli 
susceptibility, the graphene susceptibility increases with 
temperatures and number of impurities. At low temper- 
atures x(T) presents a small upturn not visible in Fig. [5] 
From the value of x(0) and i|4.3[l we obtain the critical in- 
teraction required for a ferromagnetic transition, which 



is shown in the lower left panel of Fig. [!|] Notice that 
the critical interaction strength for ferromagnetism de- 
creases as the vacancy concentration increases indicating 
that disorder favors a ferromagnetic transition. 

Using H4.4f> and the results of the previous sections we 
can also calculate the critical value for an antiferromag- 
netic transition. The result is shown in the lower right 
panel of Fig. [5J In contrast with the ferromagnetic case, 
the antiferromagnetic instability is suppressed by disor- 
der, requiring a large value of the electron-electron inter- 
action. Notice that the value of the critical ferromagnetic 
coupling is always bigger than the antiferromagnetic one, 
indicating that at half-filling the graphene lattice is more 
susceptible to antiferromagnetic correlations. This result 
is consistent with an old proposal by Linus Pauling that 
graphene should be a resonant valence bond (RVB) state 
with local singlet correlations^. 



P 0.8 
5 0.6 




,'/ 



0.2 -/ 





= 




n 


= 0.1 




n 


= 0.01 




n 


= 0.001 






0.05 0.1 0.15 0.2 



FIG. 9: (color on line) Top: Dependence of x(T) on the 
temperature and on the impurity dilution rii. Bottom: De- 
pendence of Uf and Uaf (in units of D) as function of n;. 



Hence, the Stoner criteria seems to be unable to ex- 
plain the ferromagnetic behavior observed experimen- 
tally. One might ask whether additional scattering 
mechanisms, such that provided by long-range electron- 
electron interactions, can modify the critical values of the 
couplings. The self-energy correction due to long-range 
electron-electron scattering is given in (|3.4|l and can be 
added to the Dyson equation for the Green's function 
and a new self-consistent density of states can be com- 
puted. This approach does not modify the value of Up 
which is determined by the low frequency behavior of 
the self-energy. In case of antiferromagnetism we find 
that indeed it leads to an increase on U\ F i but the re- 
sult is non-conservative since the integral over the den- 
sity of states gives a smaller value than (1 — rii). There- 
fore, we find from these calculations and previous ones^S 
that graphene is not particularly susceptible to ferromag- 
netism. 
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V. MAGNETO-TRANSPORT 

The description of the magneto-transport properties of 
electrons in a disordered honeycomb lattice is complex 
because of the interference effects associated with the 
Hofstadter problem 5 ^. As in the previous sections, we 
simplify our problem by describing the electrons in the 
honeycomb lattice as Dirac fermions in the continuum. 
A similar approach was considered by Abrikosov in the 
quantum magnetoresistance study of non-stoichiometric 
chalcogenides*&. In the case of graphene, the effective 
Hamiltonian describing Dirac fermions in a magnetic field 
(including disorder) can be written as: H = Ho+Hi with 

#0 = -V F a A- id i + e M r )} , (5-1) 

i— x,y 

where, in the Landau gauge, (A x , A y , A z ) = (— By, 0,0) 
is the vector potential for a constant magnetic field B in 
the z— direction, Oi is the i = x,y,z Pauli matrix, and 



H t = vY,Kr-rj)I. 

j=i 



(5.2) 



The formulation of the problem in second quantization 
requires the solution of Hq, which is sketched in Ap- 
pendix [BJ The field operators are defined as (see Ap- 
pendix [5] for notation; the spin index is omitted for sim- 
plicity): 



W(r) 



n.k.a 



o 

<*>(y) 



Cfc -1 



K+x(y 



-ki%) 



Ck,; 



(5.3) 



The sum over n = 0,1,2,..., is cut off at an n given 
by J3(l,7io) = D. In this representation Hq becomes 
diagonal, leading to Green's functions of the form (in 
Matsubara representation): 



Go(k, n, a; ia>) 



1 



ilu — E(a, n) 



(5.4) 



is effectively /c-independent, and E(a, —1) = is the zero 
energy Landau level. The part of Hamiltonian due to the 
impurities is written as: 



A', 



,-ixj(p-k) 



3=1 p,k 



n,a 

7i, a * 

E kltniVj ~ P l %)(t>m(yj ~ kl%) + aX^n+iiyj - pl%)(f> m +l(yj ~ M|)]4,n,a c k,m, 



n,m,a,A 



(5.5) 



r 



Equation (|5.5(l describes the elastic scattering of electrons 111 Af) : 
in Landau levels by the impurities. It is worth noting that 
this type of scattering connects Landau levels of negative 
and positive energy. 



(G(p,n,a;iu;ri, . . .r Ni )) = G(p,n,a;iu) 



= L 



-2Ni 



A. The full self-consistent Born approximation 



n / 

3=1 J 



G(p,n,a;iu;ri,...r Ni ). (5.6) 



In order to describe the effect of impurity scatter- 
ing on the magnetoresistance of graphene, the Green's 
function for Landau levels in the presence of disorder 
needs to be computed. In the context of the 2D elec- 
tron gas, an equivalent study was performed by Otha 
and Ando^iSiSiusing the averaging procedure over im- 
purity positions of Duke 57 . Here the averaging proce- 
dure over impurity positions is performed in the standard 
way, namely, having determined the Green's function for 
a given impurity configuration (ri , . . . r pf t ) , the position 
averaged Green's function is determined from (as in Sec. 



In Sec. Ill Al the averaging involved plane wave states; 
in the presence of Landau levels the average over im- 
purity positions involves the wave functions of the one- 
dimensional harmonic oscillator. In the averaging proce- 
dure we have used the following identities: 

J dy<t> n {y - pl 2 B )<j) m {y - pl%) = 5 n<m , (5.7) 
J dp<t> n {y - pl%)(t> m (y - pl%) = ■ (5-8) 
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FIG. 10: (color on line) Electronic density of states in a 
magnetic field for different dilutions and magnetic field. The 
non-disordered DOS and the position of the Landau levels 
in the absence of disorder are also shown. The two arrows in 
the upper panel show the position of the renormalized Landau 
levels (see Fig llH given by the solution of Eq. (15.1411 . The 
energy is given in units of lo c = E(l, 1). 



After a lengthy algebra, the Green's function in the pres- 
ence of vacancies, in the FSBA, can be written as: 

G(p,n,a;w + + ) = [u - E(n, a) - Ei(w)] -1 ,(5.9) 

G(p,-l;w + 0+) = [w-EaH] -1 , (5.10) 

where 

Ei(w) = -miZiuj)}- 1 , (5.11) 
E a (w) = -ni[g c G(p,-l;u + 0+)/2 + Z(w)]-\$.12) 
Z(lu) = g c G(p, -1; u; + 0+)/2 

+ . 9c ^G(p,n, a ; W + 0+)/2, (5.13) 

n,cc 

and g c = A c /(2ti1 2 b ) is the degeneracy of a Landau level 
per unit cell. One should notice that the Green's func- 
tions do not depend upon p explicitly. The self-consistent 
solution of Eqs. (£33), i|OUjl . i|QT]l . (|5~T2|) and (|5~T3|) 

gives density of states, the electron self-energy, and the 
renormalization of Landau level energy position due to 
disorder. 

The effect of dilution in the density of states of Dirac 
fermions in a magnetic field is shown in Fig. 1101 For 
reference we note that E(l, 1) = 0.14 eV, for B = 14 
T, and E(l, 1) = 0.1 eV, for B = 6 T. From Fig. 
we see that given an impurity concentration the effect 
of broadening due to vacancies is less effective as the 
magnetic field increases. It is also clear that the position 
of Landau levels is renormalized relatively to the non- 
disordered case. The renormalization of the Landau level 
position can be determined from poles of l|5.10|) : 



Of course, due to the importance of scattering at low 
energies, the solution to Eq. I|5.14|l does not represent 
exact eigenstates of system since the imaginary part of 
the self-energy is non- vanishing, however these energies 
do determine the form of the density of states, as we 
discuss below. 

In Fig. ^2 the graphical solution to Eq. (|5.14() is given 
for two different energies (E(—l, n), with n = 1,2), being 
clear that the renormalization is important for the first 
Landau level. This result is due to the increase of the 
scattering at low energies, which is present already in the 
case of zero magnetic field. The values of ui satisfying Eq. 
(|5.14|) show up in density of states as the energy values 
where the oscillations due to the Landau level quantiza- 
tion have a maximum. In Fig. 1101 the position of the 
renormalized Landau levels is shown in the upper panel 
(marked by two arrows), corresponding to the bare en- 
ergies E(—l,n), with n = 1,2. The importance of this 
renormalization decreases with the reduction of number 
of vacancies. This is clear in Fig. 1 101 where a visible shift 
toward low energies is evident when m has a small 10% 
change, from rij = 10~ 3 to n, = 9 x 10 -4 . 
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uj - E(a,ri) - ReS(o>) = 0. 



(5.14) 



FIG. 11: (color on line) Self-consistent results for Ei(w) (top) 
and Ez(w) (bottom). The energy is given in units of uj c = 
E(l, 1). In the left panels we show the intercept of ui — E(a, n) 
with ReE(tj) as required in I5.14H . 



B. Calculation of the transport properties 

The study of the magnetoresistance properties of the 
system requires the calculation of the conductivity ten- 
sor. In terms of the field operators, the current density 
operator j is given bj^: 

j = v F e[^{x 1 y)a x ^{x,y) 1 ^{x,y)a y ^{x 1 y)], (5.15) 

leading to current operator in the Landau basis written 
as: 
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J x = v F e^2~^= cl ^Ckfi^ + 4 )0)a c p _i +Uire ^[\{l-S nfi )cl nta c p ^ hX + acl n a c p ^ +1A ] , (5.16) 



C p,-l C fc,0,a C p,0,Q C P,-l 



"Fe ^ - [-A(l - d„ j0 )4 )n)a ,Cp, n _i,A + a4 inia ,Cp, n+ i !A ] -(5.17) 



As in Sec. IIIII we compute the current-current correlation function and from it the conductivity tensor is derived. 
The diagonal component of the conductivity tensor a xx (ui,T) is given by (with the spin included): 



a xx {u,T) 



4(v F e) 2 1 



dc 



2-k1 2 b uj j_ m 2tt 



\f{e + u)-f(e)] 



]- ^[ImG(-l; e + iO+)ImG(0, a; e + uj + i0 + ) 



ImG(0, a; e + iO+)ImG(-l; e + w + i0 + )] + j ^ ImG(n, a; e + iO + )lmG{n - 1, A; e + u + i0 + ) 



n>l,a,A 



- lm G(n, a;e + iO + )lmG(n + l,A;e + u> + i0 + ) 



n>0,a,\ 



(5.18) 



and the off-diagonal component a xy (uj,T) of the conductivity tensor is given by 
2(v F e) 2 1 



+ yu; + £0 + )ImG(-l; e + i0 + ) 



ReG(-l; e + juj + iO+)ImG(0, a; e + i0 + )} + 1 [ ReG K "5 e + 7^ + «0+)ImG(n - 1, A; e + i0+) 



A,n>l 



ReG(n - 1, a; e + yui + iO+)ImG(n, A; e + i0 + )}] 



(5.19) 



If we neglect the real part of the self-energy, and assume 
IxaSi(uj) = r = constant (i = 1,2), and let uj — > 0, Eq. 
(|5.18|l reduces to Eq. (85) in Ref. HU, if we further 
assume the case E(l, 1) ^> T then Eq. (|5.19|) reduces to 
Eq. (86) of the same reference. 

As in Sec. IIIII it is instructive to consider first the case 
uj,T — > 0, which leads to (a xx (0,0) = Co): 



00 



+ 



e 2 2 [ImI] 1 (0)/ImE 2 (0) - 1 
~hn [ 1 + (lm£i(0)/w c ) 2 

np + 1 

no + l + (Im£i(0)/w c ) s 



(5.20) 



When ImEi(0) ~ ImS 2 (0) and lu c > ImSi(0) (or n > 
ImEi(0)/w c ), with ui c = E(0, 1) = y/2v F /l%, Eq. (|5^U|) 
reduces to: ctq — 2/7r(e 2 //i), which is identical to the 
result l|3.11fl m the absence of the field found in Sec. 
IIIII This result was obtained previously by Ando and 
collaborators using the second order self-consistent Born 



approximation 



on. go 



However, in the FSBA it is required 



that the above conditions be satisfied for this result to 
hold. From Fig. ^|we see that the above conditions hold 
approximately over a wide ranges of field strength. 

Because the d.c. magneto-transport properties have 
been measured graphene samples^ subjected to a gate 
potential (allowing to tune the electronic density), it is 
important to compute the conductivity kernel, since this 



has direct experimental relevance. In the the case to — > 
we write the conductivity o~ xx (0,T) as: 



(X xx (0,T) 



Till 



de d -M KB (e), 
oe 



.21) 



where the conductivity kernel Kb(c) is given is Appendix 
IA1 The magnetic field dependence of kernel Kb(£) is 
shown in Fig. ^| Observe that the effect of disorder is 
the creation of a region where Kb{c) remains constant 
before it starts to increase in energy with superimposed 
oscillations coming from the Landau levels. The same 
effect, but with the absence of the oscillations, was iden- 
tified at the level of the self-consistent density of states 
plotted in Fig. [IJB Together with a xx (0,T), the Hall 
conductivity a xy (0,T) allows the calculation of the resis- 
tivity tensor: 



Pxx 
Pxy 



a 2 + a 2 

^ xx 1 xy 



a 2 + a 2 

xx 1 xy 



(5.22) 



Let us now focus on the optical conductivity, <j xx {uj). 
This quantity can be probed by reflectivity experiments 
on the sub-THz to Mid-IR frequency ranged We depict 
the behavior of Eq. (|5.18|l in Fig. ^5] for different mag- 
netic fields. It is clear that the first peak is controlled by 
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FIG. 12: (color on line) Kernel of the conductivity (in units 
of irh/e 2 ) as function of energy for different magnetic fields 
and for a dilution concentration of rii = 0.001. The horizontal 
line represents the universal limit ixhoa/e 2 — 2. 



the 25(1, 1) — 25(1, —1), and we have checked it does not 
obey any particular scaling form as function of to/B. On 
the other hand, as the effect of scattering becomes less 
important the high energy conductivity oscillations start 
obeying the scaling lj / \/B~, as we show in the lower right 
panel of Fig. 
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FIG. 13: (color on line) Optical conductivity (in units of 
nh/e 2 ) at 10 K as function of the energy for different magnetic 
fields and for a dilution concentration of m = 0.001. The ver- 
tical arrows in the upper left panel, labeled a, b, and c, repre- 
sent the renormalized E(l, 1) - E(-l, 0), E(2, 1) - E(-l, 0), 
and E(l, 1) — E(l, —1) transitions. 



VI. EXTENDED DEFECTS 



Self-doping in the absence of electron-hole 
symmetry. 



The standard description of a graphene sheet, follow- 
ing the usual treatment of the electronic band structure 
of graphitei*2iiSSiSi ; assumes that in the absence of in- 
terlayer interactions the electronic structure of graphite 
shows electron-hole symmetry. This can be justified us- 
ing a tight binding model by considering only hopping 
between ir orbitals located at nearest neighbor Carbon 
atoms. Within this approximation it can be shown that 
in certain graphene edges^iS one would find a flat sur- 
face bancfe Disclinations (a pentagonal or heptagonal 
ring) can also lead to a discrete spectrum and states at 
zero energy^? 7 . Other types of defects, like a combina- 
tion of a five-fold and seven-fold ring (a lattice disloca- 
tion) or a Stone- Wales defect (made up of two pentagons 
and two heptagons) also lead to a finite density of states 
at the Fermi level 68 i 69 i 70 . 



Band structure calculations show that the electronic 
structure of a single graphene plane is not strictly sym- 
metrical around the energy of the Dirac points^. The ab- 
sence of electron-hole symmetry shifts the energy of the 
states localized near impurities above or below the Fermi 
level, leading to a transfer of charge from/to the clean 
regions to the defects. Hence, the combination of local- 
ized defects and the lack of perfect electron-hole symme- 
try around the Dirac points leads to the possibility of 
self-doping, in addition to the usual scattering processes 
whose influence on the transport properties has been dis- 
cussed in the preceding sections. 

Point defects, like impurities and vacancies, can nu- 
cleate a few electronic states in their vicinity. Hence, a 
concentration of m impurities per Carbon atom leads to 
a change in the electronic density of the regions between 
the impurities of order m . The corresponding shift in the 
Fermi energy is ep — v-p^/nl- In addition, the impurities 

lead to a finite elastic mean free path, Z c i as ~ an i 1 ^ 2 , 
and to an elastic scattering time r c i as ~ (v-prii) , in 
agreement with the FSBA calculation in the preceding 
sections. Hence, the regions with few impurities can 
be considered low-density metals in the dirty limit, as 

T c7as - £ F- 

Extended lattice defects, like edges, grain boundaries, 
or microcracks, are likely to induce the formation of a 
number of electronic states proportional to their length, 
L/a, where a is of the order of the lattice constant. 
Hence, a distribution of extended defects of length L at 
a distance proportional to L itself gives rise to a con- 
centration of L/a carriers per unit Carbon in regions of 
order (L/a) 2 . The resulting system can be considered a 
metal with a low density of carriers, n carr ; cr oc a/I per 
unit cell, and an elastic mean free path Z c i as ~ L. Then, 
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we obtain: 



c F 



t-'F 



Tola 



(6.1) 



and, therefore, (r c i as ) _1 <C £f when a/L <C 1. Hence, 
the existence of extended defects leads to the possibil- 
ity of self-doping but maintaining most of the sample 
in the clean limit. In this regime, coherent oscillations 
of the transport properties are to be expected, although 
the observed electronic properties will correspond to a 
shifted Fermi energy with respect to the nominally neu- 
tral defect-free system. 



B. Electronic structure near extended defects. 
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We describe the effects that break electron-hole sym- 
metry near the Dirac points in terms of a finite next- 
nearest neighbor hopping between 7r orbitals, t', in 
(|2.5|) . From band structure calculations^, we expect 
that \t'/t\ < 0.2. We calculate the electronic structure 
of a ribbon of width L terminated at zigzag edges, which 
are known to lead to surface states for t' = 0. The trans- 
lational symmetry along the axis of the ribbon allows us 
to define bands in terms of the wavevector parallel to this 
axis. In Fig. ^] we show the bands closest to e = for 
a ribbon of width 200 unit cells and different values of 
t' /t. The electronic structure associated to the interior 
region (the continuum cone), projected in Fig. 1141 is not 
significantly changed by t'. The localized surface bands, 
extending from fen = (27r)/3 to ku — — (27r)/3, on the 
other hand, acquires a dispersion of order t' ( for a per- 
turbative treatment of this effect, see ref. |Z3)- Hence, 
if the Fermi energy remains unchanged at the position of 
the Dirac points (eDirac = — 3f'), this band will be filled, 
and the ribbon will no longer be charge neutral. In order 
to restore charge neutrality, the Fermi level needs to be 
shifted down (for the sign of t' chosen in the figure) by 
an amount of order t' . As a consequence, some of the 
extended states near the Dirac points are filled, leading 
to the phenomenon of self-doping. The local charge as 
function of distance to the edges, setting the Fermi en- 
ergy so that the ribbon is globally neutral. Note that the 
charge transferred to the surface states is very localized 
near the edges of the system. 



FIG. 14: Bands closest to the Dirac point of a graphene 
ribbon of 200 unit cells width. Top: t' = 0. Center: t' = 
-0.lt. Bottom: t' = -0.M 
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FIG. 15: Charge as function of the distance to the edge (in 
units of a) of a ribbon of width 300 a. Full line: t [t = —0.2. 
Broken line: t' /t = —0.1. the Fermi energy is shifted upwards 
by 0.054f for t' = -0.lt, and 0.077t for t' = -0.2i. 



Electrostatic effects. 



The charge transfer discussed in the preceding subsec- G f interaction between electrons: 
tion is suppressed by electrostatic effects, as large devi- 
ations from charge neutrality have an associated energy 

cost. In order to study these charging effects we add to Hi = \^ Ui^rnrij , (6.2) 

the free-electron Hamiltonian (|2.5|l the Coulomb energy j ,• 
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FIG. 16: Top: Self-consistent charge density (continuous 
line) and electrostatic potential (dashed line) of a graphene 
ribbon with periodic boundary conditions along the zig-zag 
edge and with a circumference of size W = 80y/3a, and a 
length of L = 960a. The parameters used are described in 
the text. The inset shows the details of the electronic density 
near the edge. Due to the presence of the edge, there is a 
displaced charge in the bulk (bottom panel) that is shown as 
a function of the width W. 



unit cell for domains of sizes 0.1 — 1/im, in agreement 
with the amount of charge observed in these systems. 
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FIG. 17: Electronic levels of a graphene ribbon with zig- 
zag edges in the presence of a magnetic field {t! = — 0.2t). 
The magnetic flux per hexagon of the honeycomb lattice is 
$ = (top), $ = 0.00025 (center), and $ = 0.0005 (bottom), 
in units of the quantum of magnetic flux, $o = h/e. The 
corresponding magnetic fields are 0T, 60T and 120T. 



where rii — X) CT ( a I o- a «>0'+&I o -&i ) a-) is the number operator 
at site Ri, and 



is the Coulomb interaction between electrons. We expect, 
on physical grounds, that an electrostatic potential builds 
up at the edges, shifting the position of the surface states, 
and reducing the charge transferred to/from them. The 
potential at the edge induced by a constant doping 6 
per Carbon atom is roughly, ~ (Se 2 /a)(W/a) (Se 2 /a is 
the Coulomb energy per Carbon), and W the width of 
the ribbon (W/a is the number of Carbons involved). 
The charge transfer is arrested when the potential shifts 
the localized states to the Fermi energy, that is, when 
t' w (e 2 /a)(W/a)5. The resulting self-doping is therefore 
5~(t'a 2 )/(e 2 W). 

We treat the Hamiltonian l|6.2l) within the Hartree ap- 
proximation (that is, we replace Hj by Hm.f. = J2i Vi n i 
where Vi = YljUi,j(nj), and solve the problem self- 
consistently for (rii)). Numerical results for graphene 
ribbons of length L — 80\/3a and different widths are 
shown in Fig.Eland Fie.nktt' /t = 0.2 and e 2 /a = 0.5t). 
The largest width studied is ~ 0.1/im, and the total num- 
ber of carbon atoms in the ribbon is s» 10 5 . Notice that 
as W increases, the self-doping decreases indicating that, 
for a perfect graphene plane (W — » oo), the self-doping 
effect disappears. For realistic parameters, we find that 
the amount of self-doping is 10~ 4 — 10~ 5 electrons per 



D. Edge and surface states in the presence of a 
magnetic field. 

We can analyze the electronic structure of a graphene 
ribbon of finite width in the presence of a magnetic field. 
The resulting tight binding equations can be considered 
as an extension of the Hofstadter problem^ to a honey- 
comb lattice with edges. The bulk electronic structure 
is characterized by the Landau level structure discussed 
in previous sections. These states are modified at the 
edges, leading to chiral edge states, as discussed in rela- 
tion to the Integer Hall Effect (IQHE) 74 . The existence 
of two Dirac points leads to two independent edge states, 
with the same chirality. In addition, Landau levels with 
positive energy should behave in an electron-like fashion, 
moving upwards in energy as their "center of gravity" ap- 
proaches the edges. Landau levels with negative energy 
should be shifted towards lower energy near the edges. 

A zig-zag edge induces also a non-chiral surface band. 
If the width of these states is much smaller than the 
magnetic field they will not be much affected by the pres- 
ence of the field. The extension of the surface states is 
comparable to the lattice spacing for most of the range 
|fc||| < (27r)/3, except near the Dirac points, so that the 
effect of realistic magnetic fields on these states is negli- 
gible. 

The finite value of the second nearest neighbor hopping 
t' modifies the Landau levels obtained from the analysis 
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of the Dirac equation. Elementary calculations (as those 
given in Appendix IB|> lead to: 



l B 4 a 2 + 2Z B V ™ , 



E±(n) = -3t' + 2l B "a ( n + - I i 

(6.4) 

with n = 0, 1, 2, 3, . . ., a = 9t'a 2 /A and 7 = Ma/2, with 
the single assumption that t 3> t' . This solution points 
out a number of interesting aspects, the most important 
of which is disappearance of the zero energy Landau level, 
made partially of holes and partially of electrons. With 
t', the electron or hole nature of the energy level becomes 
unambiguous, and half of the original zero energy Lan- 
dau level (with t' = 0)moves down in energy (relatively 
to the Fermi energy) and the other half moves up. In 
addition, the level spacing for electron and hole levels 
becomes unequal. 
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FIG. 18: (Color online) Wave function of the lowest Landau 
levels for different momenta. Black: n=0. Red: n=l. Green: 
n—2. Blue: n=3. 

The presence of a magnetic field acting on the ribbon 
does not break the translational symmetry along the di- 
rection parallel to the ribbon, that allows us to discuss 
the electronic structure in terms of the same bands cal- 
culated in the absence of the magnetic field. Results 
for the ribbon analyzed in Fig. 1141 for different magnetic 
fields are shown in Fig. El The "center of gravity" of 
the wavefunctions associated to the levels moves in the 
direction transverse to the ribbon as the momentum is 
increased. The results show the bulk Landau levels, and 
their changes as the wavefunctions approach the edges. 
The surface band is practically unchanged, except for 
small avoided crossings every time that it becomes de- 
generate with a bulk Landau level. The results show 
quite accurately the expected scaling e„ oc ±y/n for the 
eigenenergies derived from the Dirac equation, with small 
corrections due to lattice effects and a finite t' (see Ap- 
pendix[BJ) . The corresponding wavefunctions for different 
bands and momenta are shown in Fig. 1181 The Landau 



levels move rigidly towards the edges, where one also find 
surface states. 

We can compute the Hall conductivity from the num- 
ber of chiral states induced by the field at the edgesZi If 
we fix the chemical potential above the n-th level, there 
are 2 x (2n + 1) edge modes crossing the Fermi level (in- 
cluding the spin degeneracy) . Hence, the Hall conductiv- 
ity is: 



e z 4e^ 
7 xy = —2 x (2n + 1) = — (n + 1/2) . 



(6.5) 



This result should be compared with the usual IQHE in 
heterostructures, in which the factor of 1/2 is absent. 
The presence of this 1/2 factor is a direct consequence 
of the presence of the zero mode in the Dirac fermion 
problem. The existence of this anomalous IQHE was pre- 
dicted long ago in the context of high energy physicsS^ 
and more recently in the context of graphen o 12 ! 14 , but 
was observed in graphene only recently by two indepen- 
dent group a 6 i n . An incomplete IQHE, with a finite lon- 
gitudinal resistivity, was observed in HOPG graphite 77 . 



E. Fractional quantum Hall effect 

While the IQHE depends only of the cyclotron energy, 
lo Ci and therefore is a robust effect, the fractional quan- 
tum Hall effect (FQHE) is a more delicate problem since 
it is a result of electron-electron interactions. The prob- 
lem of electron-electron interactions in the presence of a 
large magnetic field in a honeycomb lattice is a complex 
problem that deserves a separate study. In this paper we 
make a few conjectures about the structure of the FQHE 
based on generic properties of Laughlin's wave functions. 

The electrons occupying the lowest Landau level are 
assumed to be in a many-body wavefunction written as 
(Ri = (xi,yi) and z = x + iy)Z& : 



* = exp(-i2m a iy j)3>(zi , . . .,z N ) . 



where 



= arctan[Im(zj — Zj)/~Re(z 



(6.6) 



/)], 



$(zi, . . . , z n) is an anti-symmetric function of the inter- 
change of two z's, and m = 0, 1, 2, The effect of the 

singular phase associated with the many-body wave func- 
tion is to introduce an effective magnetic field B* given 
by: 

B* = B — Vo(r)/e = B - 27r(2m)p(r)/e , (6.7) 
where the gauge field a(r) is given by 

(6.8) 



a(ri) = 2m'y y j V(ri)a^j , 



and p(r) is the electronic density. The procedure outlined 
above is called flux- attachment and leads to appearance 
of composite fermions. These composite particles do not 
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feel the external field B but instead an effective field B*. 
Therefore, the FQHE of electrons can be seen as an IQHE 
of these composite particles. 

Given an electronic density S, we may define an effec- 
tive filling factor p* for the composite particles as: 



P 



2nS 



(6.9) 



In the lowest Landau level the electron filling factor is 

2tt6 



P : 



eB 



and combining Eqs. I|6.9|) an d jSUHjl we obtain: 



2mp* + 1 



so, we can write: 



B* = B(l - 2mp) . 



(6.10) 



(6.11) 



(6.12) 



The crucial assumption in the case of graphene is that 
the effective p* associated with the integer quantum Hall 
effect of composite particles has the form given in Q6.5JI . 
that is (spin ignored): 



P 



(2n+l) 



n = 0, 1,2,3, . 



(6.13) 



(the effective field B* is such that the system has one 
or more filled composite particle Landau levels, and the 
chemical potential lies between two these) leading to a 
quantized Hall conductivity given by: 



2n+ 1 



2e 2 



> xy 



2m(2n+l) + l h 



(6.14) 



For n — 0, one obtains the so-called Laughlin sequence: 
Qxy — l/(2m+ l)(2e 2 /ft), and for m = we recover 
(|6.5|) . This argument shows that Jain's sequence is quite 
different from that of the 2D electrons gasS. 

As in the case of the IQHE, the FQHE can be thought 
in terms of chiral edge states, or chiral Luttinger liquids, 
that circulate at the edge of the sample 80 . One can see 
the IQHE and FQHE as direct consequences of the pres- 
ence of these edge states. Because of their chiral nature, 
edge states do not localize in the presence of disorder 
and hence the quantization of the Hall conductivity is 
robust. In graphene, as we have discussed previously, 
zig-zag edges support surface states that are non-chiral 
Luttinger liquids. We have recently shown that electron- 
electron interactions between chiral Luttinger liquids and 
non-chiral surface states can lead to instabilities of the 
chiral edge modes leading to edge reconstruction^ and 
hence to the destruction of the quantization of conduc- 
tivity. We also have shown that this edge reconstruc- 
tion depends strongly on the amount of disorder at the 
edge of the sample. While this effect is not strong in the 
IQHE (because the cyclotron energy is very large when 
compared with the other energy scales), it makes the ex- 
perimental observation of the FQHE in graphene very 
difficult. 



VII. CONCLUDING REMARKS 

To summarize, we have analyzed the influence of local 
and extended lattice defects in the electronic properties 
of single graphene layer. Our results show that: (1) Point 
defects, such as vacancies, lead to an enhancement of the 
density of states at low energies and to a finite density 
of states at the Dirac point (in contrast to the clean case 
where the density of states vanishes); (2) Vacancies have 
a strong effect in the Dirac fermion self-energy leading 
to a very short quasi-particle lifetime at low energies; (3) 
The interplay between local defects and electron-electron 
interaction lead to the existence of a minimum in the 
imaginary part of the electron self-energy (a result that 
can be measured in ARPES); (4) The low temperature 
d.c. conductivity is a universal number, independent on 
the disorder concentration and magnetic field; (5) The 
d.c. conductivity, as in the case of a semiconductor, in- 
creases with temperature and chemical potential (a result 
that can be observed by applying a bias voltage to the 
system); (6) The a.c. conductivity increases with fre- 
quency at low frequency and at very low impurity con- 
centrations can be fitted by a Drude-like model; (7) The 
magnetic susceptibility of graphene increases with tem- 
perature (it is not Pauli-like, as in an ordinary metal) and 
is sensitive to the amount of disorder in the system (it 
increases with disorder); (8) Within the Stoner criteria 
for magnetic instabilities we find that graphene is very 
stable against magnetic ordering and that the phase dia- 
gram of the system is dominated by paramagnetism; (9) 
In the presence of a magnetic field and disorder, the elec- 
tronic density of states shows oscillations due to the pres- 
ence of Landau levels which are shifted from their posi- 
tions because of disorder; (10) The magneto-conductivity 
presents oscillations in the presence of fields and that 
their dependence with chemical potential and frequency 
are rather non-trivial, showing transitions between differ- 
ent Landau levels; (11) Extended defects, such as edges, 
lead to the effect of self-doping where charge is trans- 
fcrcd from/to the defects to the bulk in the absence of 
particle-hole symmetry; (12) The effect of extended de- 
fects on transport is very weak and that electron scatter- 
ing is dominated by local defects such as vacancies; (13) 
The quantization of the Hall conductance in the IQHE is 
anomalous relative to the case of the 2D electron gas with 
an extra factor of 1/2 due to the presence of a zero mode 
in the Dirac fermion dispersion; (14) We conjecture that 
the FQHE in graphene has a sequence of states which is 
very different from the sequence found in the 2D electron 
gas and we propose a formula for that sequence. 

The results and experimental predictions made in this 
work are based on a careful analysis of the problem 
of Dirac fermions in the presence of disorder, electron- 
electron interactions and external fields. We use well 
established theoretical techniques and find results that 
agree quite well with a series of amazing new experi- 
ments in graphene^SiL&ifliii. The main lesson of our 
work is that graphene presents a completely new electro- 
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dynamics when compared to ordinary metals which are 
described quite well within Landau's Fermi liquid the- 
ory. In this work, we focus on the effects of disorder and 
electron-electron interaction and have shown that Dirac 
fcrmion respond to these perturbations in a way which is 
quite different from ordinary electrons. In fact, graphcnc 
is a non-Fermi liquid material where there is no concept 
of an effective mass and, therefore, a system where Fermi 
liquid concepts are not directly applicable. A new phe- 
nomenology, beyond Fermi liquid theory, has to be de- 
veloped for this system. Our work can be considered a 
first step in that direction. 
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APPENDIX A: 9(w,e), AND ifs(e) 
In the calculation of a(u,T) and a(0,T) we defined: 



6(u;,e) = 



Sl=±l,S2 
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log 



=±i 
(D - Sl Af - 
A 2 + B 2 



arctan 
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where and cosa = (C 2 - E 2 )j(C 2 + E 2 ). 

M = (C 2 + E 2 - A 2 - B 2 )/V, V = —M , 

N = 2( Si A- s 2 C){A 2 +B 2 )/V, 

P = -2( Sl A - s 2 C)(A 2 +B 2 )/V, 

V = (A 2 + B 2 -C 2 - E 2 ) 2 

+ 4(A 2 + B 2 )(C 2 - s x s 2 AC) , 

+ 4(C 2 + E 2 )(A 2 - Sl s 2 AC) , 

A = e + w-Re£(e + w), 

B = Im£(e + w), 

C e ReE(e) , ^ ^ e ma g ne to-transport properties, a xx (0,T) given 

E = ImE(e) , by gq J5 21|| . depends on the kernel Kb{x), which is 

(A2) defined as: 



K B (x) 
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A,n>l 
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Iml^a:) 
ReS 2 (x)] 2 + [ImS 2 (x) 

ImSi(a;) 



ImKi(i) 

E(a,0) -ReEi(x)] 2 + [ImSi(a:)] 2 
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[x - E(a, n) - ReSi(x)] 2 + [ImSi(a;)] 2 [x - E(X, n-l)- ReZ^x) 



[Im£i(x)p 



(A3) 



APPENDIX B: THE DIRAC EQUATION IN A with L the size of the system in the x direction. After 
MAGNETIC FIELD straightforward manipulations, the eigenproblem reduces 

The Hamiltonian l|5.1|l can be solved using a trial 
spinor of the form: 
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to: 



where 



is an eigenfunction of (|B2|I with eigenvalue E(a,n) = 
+ 1, with a = ±1, and <f> n ( n = 0, 1, 2, . . .) 
(B2) the n eigenfunction of the usual ID harmonic oscillator. 

In addition, there exists a zero energy mode whose eigen- 
function is given by: 



V2h 



V2h 



■(y + i%d v ) 
-{y-ildy) 



(B3) 
(B4) 



with the magnetic length defined as l 2 B — l/(eB). For 
the case of E ^ 0, it is simple to see that the spinor 



±( My) 
y/2 V a <f>n+i(y) 



o 

i>o(y) 



(B6) 



that completes the solution of the original eigenproblem. 
As in the more conventional Landau level problem, the 



degeneracy of each level is L B/4>o, with 



h/e the 



(B5) quantum of flux. 
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